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We report the first observation of B° — > J/ipil and evidence for B^ —* J/ipr)' . These re- 
sults are obtained from 23.6 fb~^ of data collected at the T(55) resonance with the Belle de- 
tector at the KEKB e'^e~ collider. We measure the branching fractions B{B^ — > J/iprj) = 
(3.32 ± 0.87(stat.)tg;^|(syst.) ± 0.42(/s)) x 10"^ with a significance of 7.3a, and B{B° ~* Jf-ipv') = 
(3.1 ± 1.2(stat.)l^ j(syst.) ± 0.38(/s)) x 10"'' with a significance of 3.8cr. 

PACS numbers: 13.25.Hw, 14.40.Nd 

The decays J/ip'q^'^ are dominated by the h used to measure directly the i?^ width difference AFg [l[ . 

CCS process, as shown in Figure [T] The J/iprj'''^ final Using SU(3) flavor symmetry, an estimate of the — > 
states are CP eigenstates, whose time distribution can be J/ip-q''-'^ branching fractions can be obtained relative to 
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FIG. 1: Dominant diagram for the processes — > J/i^V 



the decay -> J/tPK° 
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where p* is the momentum of J/^ or 77^'' in the 
rest frame of the or _Bj]. Here, i/ip ~ 37° is 
the quark mixing angle in the flavor basis with j] = 
-^{uu + dd) cos (j)p — ss sin (f>p. We then expect B{B'^ 

J/V'ryC)) - 3.45(4.88) x 10-^ using the value B{B°i 
J/^PK^) = 8.71 X 10-4 [4]. To date there is only the up- 
per hmit B{B° J/^rf) < 3.8 x lO'^ at 90% confidence 
level [55. Recently, a large sample of B^ mesons was 
produced from e+e" collisions at the T(5S') resonance, 
where final states with photons can be reconstructed with 
low background. 

In this letter, we report measurements of fully re- 
constructed i?° J/i^V and B^ J/i^v' decays us- 
ing a (23.6 ± 0.3) fb~^ data sample collected with the 
Belle detector at the KEKB asymmetric-energy e+e^ 
collider operated at the T(5iS') resonance. The 
beam energy in the center-of-mass (CM) frame i?bcam 
is measured to be 5433.5 ± 0.5 MeV using T(5S') 
T(lS')7r+7r-, T{1S) tj,- decays Q. The total hh 

cross section at the T(5S') energy has been measured 
to be crj^'^^^ = (0.302 ± 0.014) nb 9]. B° events are 

produced from _Bs*^i3s*'' pairs where three B^ produc- 
tion channels B°B'^, B*B'^, B*B* are kinematically 
allowed. The total fraction of Bl ' Bl ' pairs in hh 



events has been measured to he fs — N 



bb 



m- Thus, the 23.6 ib"^ data sample con- 
tains a total of 2.78 million mesons. The frac- 
tions of -B^ production channels are measured to be 
fe^S. = Nb.s,/N^msm = (90.ll^:«±0.2)%, /^.^o = 
Ns^so/N^M^M = (7.3+S5:33i^±0.1)% 0. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 



ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect mesons and identify muons (KLM). The 
detector is described in detail elsewhere 

Charged tracks are required to originate within 0.5 cm 
in the radial direction and within 5 cm in the beam di- 
rection, with respect to the interaction point. Electron 
candidates are identified by combining information from 
the ECL, the CDC (dE/dx), and the ACC. Muon candi- 
dates are identified through track penetration depth and 
hit patterns in the KLM system. Identification of pions is 
based on combining information from the CDC {dE/dx), 
the TOF and the ACC. 

Two oppositely charged leptons {I — e or fi) and 
bremsstrahlung photons lying within 50 mrad of or 
tracks are combined to form a J /ip meson. The leptons 
are required to be positively identified as electrons or 
muons. The invariant mass is required to lie in the ranges 
-0.150 GeV/c2 < A4e(^) - mj/^ < 0.036 GeV/c^ and 



-0.060 GeV/c2 < Af^ 



ij/^ < 0.036 GeV/c^ where 



mj/^ denotes the nominal J /if: mass, and A/ge(-j,) and 
A/pp are the reconstructed invariant masses for 6+6^(7) 
and /i^/i^, respectively. 

Photon candidates are selected from ECL showers not 
associated with charged tracks. An energy deposition 
with a photon-like shape and an energy greater than 50 
MeV is required. 7r° candidates are selected by combin- 
ing two photon candidates with an invariant mass in the 
range 115 MeV/c^ < M^^ < 155 MeV/c^ . 

Candidate rj mesons are reconstructed in the 77 
and iT^Tr~Tr'^ final states. We require the invariant 
mass to be in the range 500 MeV/c^ < M~^j < 
575 MeV/c2 ([-3.5(j, 2.0cr]) and 535 MeV/c^ < 
M^+^-^o < 560 MeV/c^ (±2.5f7). In the 77 final state, 
we require that |cos6'doc| < 0.9, to reduce background 
from combinatorial 7's, where 0dec is the angle between 
the 7 and 77 lab momentum direction in the rj rest frame. 

Candidate rj' mesons are reconstructed in the riTr~^Tr~ 
and /9°7 channels. The rj candidates are selected in the 
same two channels as above. Thus, there are three sub- 
channels for 77' reconstruction, p^'s are selected from 
oppositely charged pion pairs satisfying 550 MeV/c^ < 
A/^+^- < 900 MeV/c^ and a heficity requirement 
Icos^hcil < 0.85 since the p° is longitudinally polarized. 
Here 6*1101 is the helicity angle of p", calculated as the 
angle between the direction of the 7r+ and the direc- 
tion opposite to the ij' momentum in the p'^ rest frame. 
We require the reconstructed rj' invariant mass to satisfy 
940 MeV/c2 < Af^, < 975 MeV/c^ (±3cr). 

We combine J/i/; and r]^'^ candidates to form B^ 
mesons. Signal candidates are identified by two kine- 
matic variables computed in the T(5S') frame: the energy 
difference AE = Eg — i?bcam and the beam-energy con- 
strained mass ATbc = \/(-Ebcam)^ — (Pij)^i where Eg and 
p*g are the energy and momentum of the reconstructed 
candidate. To improve the AE and A/bc resolutions, 
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mass-constrained kinematic fits are applied to J/'0, r?*-'-* 
and TT^ candidates. We retain meson candidates with 
|A£:| < 0.4 GeV and Mtc > 5.25 GeV/c^ for further 
analysis. If there are multiple candidates in a single 
event, we choose the candidate that minimizes the sum 
of the x^'s of the mass-constrained fits. 

To suppress the two-jet-like continuum background 
from e~^e~ ^ qq {q = u,d,s,c), we require the ratio 
of second to zeroth Fox- Wolfram moments IjJ to be less 
than 0.4. This requirement is optimized by maximizing a 
figure-of- merit Ns/\/Ns + Nb, where Ns is the expected 
number of signal events and Nb is the number of con- 
tinuum events, in the B*B* signal region. The contin- 
uum background is modeled by a second-order polyno- 
mial in AE and an ARGUS function [l3| in Mbc. We 
study the continuum background in a J/V' sideband de- 
fined as 2.5 GeV/c^ < Mj,^ < 3.4 GeV/c^, excluding 



the region -0.200(-0.080) GeV/c^ < Mj/^ - mj/^ < 
0.048 GeV/c^ for J/ijj e+e~{ti+^i-) in data. The 
shapes for continuum backgrounds for each sub-channel 
are determined from fits to J/ip sideband data with a 
relaxed lepton identification requirement. The yields for 
continuum background are determined from the yields 
in J/ip sideband data with all selections other than the 
J/i/j mass cut as in the nominal selection. Scale factors 
are determined from the ratio of the J/tp selection area 
to the J/ip sideband area by fitting Monte Carlo (MC) 
J/^ mass distributions. For continuum J/tp backgrounds 
that may not be well modeled in MC, we studied the off- 
resonance data and obtained the fractions of real J/-0's 
in continuum, which are applied as a correction to the 
continuum background yields. 

The remaining background is from BB [B = 
Bg, B^) events with one B meson decays to a 

final state with J/'ip (denoted J/ipX). We use a MC 
sample generated at the T(5S') resonance that includes 
all known B — + J/ipX processes to estimate this back- 
ground, and find that the dominant contribution is from 
Bjj, B^ — * J/ip + strange mesons. This background does 
not peak in either AE and Mbc and is described by 
an exponential function in AE and an ARGUS func- 
tion in Mbc, with shapes determined from MC. For 
the channel J / tprj' {p'^ , the size of the J/tpX back- 
ground is comparable to the expected signal yield in the 
B*Bl signal region, so we use a rj' sideband defined as 
0.90 GeV/c2 < M,,- < 0.93 GeV/c^ and 0.99 GeV/c^ < 
Mn' < 1.02 GeV/c^ to study it. We fit the -q' side- 
band data to obtain the J/ipX background shape for the 
J/iprj' (p'^j) channel, where the continuum contribution 
is fixed from the J/ip sideband. 

We categorize events into two rj and three 77' sub- 
channels. For the J/tpV and J/ipv' niodes, we per- 
form a simultaneous unbinned maximum likelihood fit 
to the AE — Mbc distributions for each group of two ij 
or three rj' sub-channels, with the signal branching frac- 
tion as a common parameter while treating the signal 
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FIG. 2: Distributions for the J/ipi] channels, (a) Mbc dis- 
tribution for the B*B* signal region with Ai5 € [—97, —9] 
MeV, and (b) A_E distribution for the B*B* signal region 
with i\fbc G [5.405,5.431] GeV/c^ The full histogram shows 
projections of fit results. The sum of all backgrounds is repre- 
sented by the blue dotted curves with the green dashed curve 
corresponding to continuum background. 
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FIG. 3: Fit projections for the combined clean J/'tprj' {tj-k'^tt') 
channels (a,b) and J/ipr/' (p'^j) channel (c,d). The projections 
are shown in the B^Bg signal region with AE £ [—98, —17] 
MeV (a,c), and with Mbc G [5.403,5.428] GeV/c^ (b,d). The 
sum of all backgrounds is represented by the blue dotted 
curves with the green dashed curve corresponding to contin- 
uum background. 



and background shapes separately for the different r]^'^ 
sub-channels. 

In the fit, the signal normalization 
production channel is parameterized as 



N 



for each 5° 
Nsig = 2 X 
We use three 

fractions fg, ; /b* - and f^o so = 1 - /g. 5. - Z^. 50 . 
In the J/tpr]' mode, because of low statistics, we only 
include the B*B* channel in the fit. The index i de- 
notes each 77^'-' sub-channel. The product Bi = B^J/tp 
l^l^)Bi{ri'^'^) is the total branching fraction to final states 
with a J/-0 and an ri'^'\ and ei is the MC reconstruction 
efficiency. The values of the weighted efficiencies Biei 
are listed in Table [J The signal shapes are from signal 
MC histograms while the means and widths of the dis- 
tributions are corrected using a B^ — > J/ipK*'^{K*~^ — > 
-ftT+Tr'^) control sample from T(4S'). The continuum back- 
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TABLE I: A summary of efficiencies, significances and branching fractions (B) for each mode. 



mode 




Significance 




B 


B° — > T/ihn('y'y) 

B° J/V'?7(7r+7r-7r") 

Total B° J/tpri 


1.34% 
0.45% 
1.79% 


7.3(7 


(3.32 ± 0.87(stat. 


)lo;28(syst.)±0.42(/,))xl0-^ 


B° ^ j/Hin^^) 
B° J/Hip'y) 

Total ^ J/V'r;' 


0.52% 
0.88% 
1.40% 


3.8cr 


(3.1 ± 1.2(stat.; 


llJ!j(syst.)±0.38(/s)) X 10-* 



TABLE IL Relative systematic errors (in %) for B{J/tpri'^'^). 



Source 


s(J/V"?) 




Signal shape calibration 


+5.8,- 


2.9 


+11.7,- 


16.8 


Beam energy 


+1.6,- 


0.0 


+4.8, - 


4.3 


MC signal shape 


+1.0,- 


2.0 


+2.6,- 


4.0 




+0.7, - 


1.5 


+4.6, - 


4.0 


Background parameters 


+0.9,- 


0.8 


+6.0,- 


5.5 


Track reconstruction 


2.5 




4.2 




Lepton identification 


4.2 




4.1 




Pion identification 


0.4 




2.3 




77(71'^) — > 77 selection 


4.1 




2.8 




B{J/tP 11) 


0.72 




0.72 




B{ti'-'^ final states) 


0.49 




2.3 




Luminosity 


1.3 




1.3 






4.6 




4.6 




/. 


+13.4,- 


13.3 


+13.4,- 


13.3 


Total 


+16.8,- 


16.1 


+21.9,- 


24.8 



ground's shapes and yields and the J/ipX background's 
shapes are fixed. The floating parameters for each fit are 
the branching fraction B{B'^ J/iprj'^'^) and the J/ijjX 
background yields for each sub-channel. 

The projections of the fit in the B*B* signal region 
are shown in Figures [2] and [3l The signal efficiencies, 
branching fractions, and significances including system- 
atic uncertainties are listed in Table HI We calculate a 
total of 14.9 ± 4.1 B° J/t/jj] events_and 10.7 ± 4.6 
B° J/ip-q' events in the T(5S') ^ B*Bl channel. 

The B^ J/ipT] decay is observed for the first time 
and evidence for the i?° — > J/ipjf decay is found. The 
significance is defined by S* = •\/21n(£max/'Co), where 
'Cmax('Co) is the likelihood value at the maximum (with 
the signal branching fraction set to zero). The signif- 
icance including systematic uncertainty is taken as the 
smallest value for each systematic variation described be- 
low. 

The muon and pion identification efficiencies from MC 
are calibrated using the J/ip — > l^l^ and D^tt^ 
control samples in data, respectively. The total system- 
atic error due to lepton identification is weighted to be 



4.2(4.1)% for J/iprj'^'^ modes. The systematic error due 
to pion identification is 0.4(2.3)% for the J/f/'Ty'^'^ modes. 

The systematic errors due to the signal shape mean and 
width corrections and background parameters are deter- 
mined by varying each parameter by its error, repeating 
the fit, and summing the shifts in branching fraction in 
quadrature. The beam energy has an error of ±0.5 MeV, 
whose systematic effect is evaluated by varying the mean 
value of Ai? and Mbc for the signal shapes simultaneously 
according to the uncertainty in the beam energy. All the 
systematic errors are summarized in Table |lll The large 
systematic errors due to fs are quoted separately in the 
final results. 

The ratio of the two branching fractions R — B{B^ 
J / ip-q') / B{B1 J/il'v) is also calculated. The statisti- 
cal errors of the two modes are combined. The common 
systematic errors due to luminosity, cross-section and fs 
cancel. Correlated systematic errors due to calibration, 
beam energy, track reconstruction and particle identifica- 
tion are treated properly by varying the numerator and 
denominator simultaneously. Other systematic sources 
in the two branching fractions are treated independently. 

In summary, we observe i?^ — > J/iprj decay with a sig- 
nificance of 7.3(7 and find evidence for B^ J/il'v' with 
a significance of 3.8(t. We measure the branching frac- 
tions B{B° J/iPt]) = (3.32±0.87(stat.)l°-g(syst.) ± 
0A2{fs)) X 10-4 and ^(5° ^ J/H) = (3-1 ± 
1.2(stat.)t[] j(syst.) ± 0.38{fs)) x lO^"*. The ratio of two 
branching fractions is measured as i? = B(B^^jjt>r)) ~ 
0.924t°;^4(stat.);j5;2o(syst.). The results are consistent 
with SU(3) expectations using the measured value of 
6(5°^ J/^ifO) ii. 
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